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Ci . Abstract 

Axisymmetric accretion disks in vicinity of a central compact body are 
Q ■ studied. For the simple models such as vertically isothermal disks as well 

5— l ' as adiabatic ones the exact solutions to the steady-state MHD (magneto- 

C/3 , hydrodynamic) system were found under the assumption that the radial 

Cw ■ components of velocity and magnetic field are negligible. On the basis of 

the exact solution one may conclude that vertically isothermal disks will 

^sj . be totally isothermal. The exact solution for the case of adiabatic disk cor- 

£> roborates the view that thin disk accretion must be highly nonadiabatic. 

Q\ _ An intermediate approach, that is between the above-listed two, for the 

l/~) ■ modeling of thin accretion disks is developed. In the case of non-magnetic 

disk, this approach enables to prove, with ease, that all solutions for the 
midplane circular velocity are unstable provided the disk is non-viscous. 
■^1- ■ Hence, this approach enables to demonstrate that the pure hydrodynamic 

<T^) turbulence in accretion disks is possible. It is interesting that a turbulent 

£/") . magnetic disk tends to be Keplerian. This can easily be shown by assum- 

ing that the turbulent gas tends to flow with minimal losses, i.e. to have 
the Euler number as small as possible. 
• y—t - 

H '. 1 Introduction 

c3 



We will consider the dynamics of axisymmetric accretion disk around a compact 
object. A successful theory of the process in question is mainly developed (see, 
e -g- 0, [2], [22], [20], [21]). Extensive use is made of simple models such as 
vertically isothermal disks as well as adiabatic ones. To estimate possible errors 
and limits associated to these models, the exact solutions to the steady-state 
MHD (magnetohydrodynamic) system were found (see Sec. [5]and Sec. [3]) under 
the assumption that the radial components of velocity and magnetic field are 
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negligible. We will also consider an intermediate case that is between the above- 
listed two axisymmetric flows (see Sec. [4]). Such an approach for the modeling 
of thin accretion disks turns out to be more flexible and efficient. In particular, 
the question of pure hydrodynamic turbulence is still an open question [7] . The 
possibility for finite disturbances to develop turbulence in the nonlinear regime 
was demonstrated by O. A. Kuznetsov in [8] and in doing so he has disproved 
the well-known arguments that pure hydrodynamic turbulence cannot be a self- 
sustaining source of viscosity in accretion disks (see, e.g., [7] and references 
therein). Using the approach of Sec. 2] we will also demonstrate that the pure 
hydrodynamic turbulence in accretion disks is possible. 

The input system of magnetohydrodynamic (MHD) equations is the follow- 
ing (see, e.g. [9j p. 126], [H], [15], [16], [17], [H], EH): 
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-pv-V$ - P Q - V • q, (3) 



<9B 

~dt 



= V x (v x B) , 
V-B = 0, 
E = -~(vxB), 



(4) 
(5) 

(6) 



where p, v, P, B, E, t, q, and $ denote the density, velocity, pressure, mag- 
netic induction field, electric field, stress tensor, heat current, and gravitational 
potential, respectively, B — |B|, v = |v|, E = pe p + 0.5pv 2 + B 2 / {Sir) denotes 
the total energy per unit volume with e p being the internal energy per unit 
mass for the plasma, Q denotes the local cooling rate [51 p. 142]. It is, mainly, 
assumed that 

$ = -G J^-, , G, M = const. (7) 

Vr z + z z 

The heat conductive flux, q, may be expressed as 

q = -A T V ■ T, 



(8) 



where At denotes the thermal conductivity, T denotes the temperature. The 
stress tensor, t, is the sum of two tensors, r = t v + T tl namely, the viscous, 
t„, and the turbulent, Tt, stress tensors: 



u v [Vv + (Vv)*] - -p,V • vl, 



0) 



r t « Mt [Vv + (Vv)*] - | (m*V • w+pK) I, (10) 

where ( ) denotes a conjugate tensor, \x v denotes the dynamic viscosity, fx t 
and k denote the turbulent viscosity and the kinetic energy of turbulence pQ, 
respectively. We will also use the viscosity n = fi v + fi t . Obviously, if the flow 
is laminar, then k — and fj, is the dynamic viscosity. 

We will not take into consideration the energy equation, since the following 
three axisymmetric flows in cylindrical coordinates, (r, 93, z), will be considered. 

1) A vertically isothermal disk, where the temperature is a pre-assigned 
value, 

dT 
P = pRT, — = 0, R = const. (11) 

2) An adiabatic disk, i.e. 

P = Kp 1 1 1,K = const. (12) 

3) An intermediate case that is between the above-listed two, in some mea- 
sure opposite, axisymmetric flows. 

Let us note that the number densities of ions and electrons at any point 
are approximately equal, and, hence, a plasma must always be close to charge 
neutrality. Even a small charge imbalance would create huge electric fields 
which would move the plasma particles so as to restore neutrality very quickly 
[5], [7]. The plasma maintains charge neutrality to a high degree of accuracy. 
However, local charge imbalances may be produced by thermal fluctuations [9]. 
To estimate their size, the Debye length is usually used [5], [9], [7j. The length 
scale of plasma dynamics must be much larger than the Debye length. For 
example, the mean density of gas in the Milky Way is a million per cubic metre 
[S]. Then, assuming that the gas temperature is close to 100° K, we find that 
the Debye length will be less than 0.7 to. We will consider the flow at the 
periphery of accretion disk and, hence, the length scale of plasma dynamics 
will be much larger than the Debye length. If, however, inside the disk there 
is a charge density, then it gives rise to an electric field outside the disc which 
is available to pull charges out of the disc [7j. Hence, in the case of at least 
steady-state flow, we may write that 

V-E = 0^ V-(vxB) = 0. (13) 

Let us introduce the following characteristic quantities: £*, /*, p„, i>*, p*, 
T*, //„, k*, and B* for, respectively, time, length, density, velocity, pressure, 
temperature, viscosity, kinetic energy of turbulence, and magnetic field. The 
following notation will also be used: 



h w _ P* „ _ ^P* p _ v%l* R _ P*vJ* 9 _ 2k» 
v*t* P*vi Bi GM IX* 3vi 



where Sh, E u , F r , and R e denote, respectively, Strouhal, Euler, Froude, and 
Reynolds numbers. For axisymmetrical flow, we have, in view of (|14p. the 



following non-dimensional system: 



dp ld(rpv r ) d(pv z ) 

<->/» — H E 1 E = U i 



dt r dr 



dz 
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dpv r 19 / 2 £„ 2 



(9 



pV r V z 



f3 r r 

J_(d_ ' 







dz 
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p d<§> 



2 dvr_ _ 2 /l d(ny) dv z 
dr 3 \r dr dz 



d ( dv r dv z \ 2/i / dv r v r 
dz \ dz dr J r \ dr r 
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dpv v 19 / E u 

Sh -bT + rd-r r { PV ^ Vr -T ¥ ' r 



d ( E u \ pv v v r E u B V B 

Tz [ pv * Vz T v 
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d fv t 
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r 
d ( dv t 



dz V dz 



/3 r 
d fVu. 



2p 



dr 



dpv z 19 / E u 

bh^r— + ~T-r pv z v r ^-B Z B 7 

dt r dr \ p 



dz 



pv 



2 f^L p2 



'dz [ EuP+ I 2 ) F r dz' 



j_ (d_ 

R'eX'dr' 



dv r dv z 



dz 
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d_ 

dz 
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P_ I ' dvj_ chh 
r V dz dr 



dv z _ 

dz 

-0k, 



1 d(rv r ) 
3 \r dr 

d 



dv z 

dz 



d. 



S>, 



dB r d (v z B r — v r B z 



-pK, 



0, 



s h 



dt dz 

dB v d{v r B v -v v B r ) d{v z B v -v v B z ) 



dt 



dr 



dz 



= 0, 



S h 



dB z 



1 dr (v r B z - v z B r ) 



dt r dr 

1 drB r dB z 
r dr dz 



= 0. 



0. 



It is, mainly, assumed that 



(16) 



(17) 



(18) 
(19) 
(20) 
(21) 
(22) 



$ = — 



\Jr 2 + z 2 



(23) 



We will consider, in general, accretion disks. Hence, it is assumed that the 
values p, iv, v Vl P, $ are even functions of z, whereas v z is an odd one. In 
such a case, in view of (|15[) - (j2"5)) . there exist two possibilities: 1) B z is an even 
function of z, whereas B v and B r are odd ones; 2) the values B v and B r are even 
functions of z, whereas B z is an odd one. We will consider the first possibility. 

If v r = v z = 0, then there exists the third possibility, namely, the values 
B v and B z are even functions of z, whereas B r is an odd one. In such a case 
the magnetic field will be unstable provided B v ^ 0. Actually, if v r ^ 0, then 
B v will be an odd function of z, since we consider the case when B z is an even 
function of z. It is very important to note that the solution such that B v (^ 0) 
is an even functions of z can not be obtained as a limiting case (namely, as 
v r — > 0) of the motion under v r ^ 0. Hence, any solution for B v (^ 0) such 
that it is not an odd function of z may be seen as unstable, as any infinitesimal 
variation, Sv r (7^ 0), gives rise to a finite response in the magnetic field. 

In the case of steady-state flow, we mainly assume that the flow is charge- 
neutral HH, i.e. 



dr (v v B z - v z B v ) d (v r B v - v v B r 



0. (24) 



rdr dz 

2 Vertically isothermal disk 

In this section we construct steady-state solutions for the system (|T5|)-([23ll pro- 
vided that 

dT 
v r = 0, — = 0. (25) 

oz 

Let us note that the equality v r — implies v z — 0. It could be easily seen from 
the steady-state version of (fT5J) . We take p* = p^RT* and, hence, we obtain 
from dH]) that 

P = pT. (26) 

The temperature in (J26I) is assumed to be a preassigned function of r. 

We intend to find an exact steady-state solution to the MHD system (Un- 
provided R e — ¥ 00, dke = 0, and 



B r = 0. (27) 

In such a case, the MHD system is reduced to the following. 

E u Bl pv%_ d( E U B*\ pd$ 

p r r ~ dr\ u+ (3 2 ) F r dr ' {M) 

d(B v B z ) 



dz 



0, (29) 



8 (E U A _ d (^ D , E U B 2 \ p 9$ 



M^J = ^r + 7Tj-^ (3o) 



d(v v B 2 ) 

dz 

dB z 



= 0, 



dz 



= 0. 



In view of (|32[). (jBTl) . and ([29]), we obtain: 

S 2 = B 2 (r) , v v = v v (r) , B v = B v (r) 
Then we obtain, instead of I 



£„ £2 pw 2 



(9 



^-t-'^ = --\E u P 
p r r or 



E u B\ + B 2 \ P 9$ 
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F r dr' 



o-i^-m 



C P = P\z=0 ■ 



Let 

By virtue of (|26| . we find from (|35|) that 



$1 



z=0 ' 



p = C„exp. 
If dM]) is valid, then 

p = C p exp 
Let, in general, 



i 



f 



i 



TF r E u Wr 2 + z 2 r 
9$ 



C p = C p (r) . 



dz 



*0. 



(31) 
(32) 

(33) 

(34) 

(35) 
(36) 



C p = C p (r), = ^.(r) ^ <&L =0 . (37) 



(38) 



(39) 



Eq. (|3"4"|) must be valid under all values of z > 0. After differentiation (|3"4"|) over 
z, in view of (|35[) we obtain 



v% p _ dp p d<f> 

r TE U ~ dr + TF r E u dr ' 



Since, in view of (|37l) . 

dp 

dr 
we find from (|40p that 



dc p a ( t-s 

dr p dr \TF r E u 



exp 



\TF r E u 



b c p {(t>- $) dr 

~TE~ u ~lfr^TF~E~ u dr~ T 2 F r E u dr' 



V ip *~>p _ dCyp Cp 



(40) 



(41) 



(42) 



After differentiation (|4"2"]l with respect to z, we obtain: 

-3^**1=0 (43) 

T 2 F r E u dz dr { ' 

Since C p 7^ 0, we obtain from (j43|) . in view of (|39|) . that 

T = const, (44) 

and, hence, 

C a = const exp / — -^— dr — — . (45) 

p P ^7 rT^ u TF r E u ) V ' 

If ((231) is valid, then </> = —\/r. Let the disk be Keplerian, i.e. 

1 



\Jrb r 
then 



(46) 



dr 1 



C p = constexp \J -^^- + ^r^J = const (47) 

Thus, the assumption that the motion is Keplerian leads to a constant density 
at the midplane. Let us consider a vortex motion, i.e. 

C 
v v = — — , C v = const, (48) 

and let, for the sake of simplicity, C v = \/ \[F~ r , then, in general, we have 



*£ < 0. (49) 

The density at the midplane, in view of (14"5"|) . will be the following. 

C n = const exp = I . (50) 

9 V \rTF r E u 2r 2 TF r E u ) V ; 

Substituting (J37)) . (|44l) . and (|45)) into ([34j) , we obtain the following equation in 
B v and B z . 

Let us assume that the flow is electrically neutral [21] to close Eq. (|5T|) . 
Hence 

^%^1 = 0. (52) 

ror 

If we have a vortex flow, i.e. (14"51) is valid, then 



B z = const. (53) 



By virtue of (f55|). we find from (j5Tj) : 



B v = — -, Cb v = const. 



(54) 



Let us note that the solution, (|54|) . for i3 v is unstable, as this solution is an 
even function of z (see Sec. [T|), and, hence, B v = 0. 

Conversely, we can see from (155)) that B^ (^ 0) is an even function of z, and, 
hence, B v is unstable (see Sec. [J). Thus, the stable solution is: B v = 0, and, by 
virtue of (|5"TT) . we find that (|53l) is valid. Assuming that the flow is electrically 
neutral [3T], we find from (|52l) that (|4"8"]l is valid, i.e. we have the vortex flow. 

The semi-thickness, H, of disk is often (e.g. [6]) defined as 



H 



C n 



pdz 7 C p = p\ z= 



z=0 ■ 



(55) 



Notice, using the exact solution ([38]) in (1551) wc find that i? — > oo provided that 
TF r E u 7^ 0. Thus, even if the value of E u be small but finite, the disk cannot 
be thin in terms of ([55]) . Instead of the exact solution, (f38| , it can be used the 
following approximation (Cf. [18], [21]) for small values of z. 



C p exp I - 



2TF r E u r 3 



— C p exp 



\{z/Hf 



H = r^rTF r E u . (56) 



Let us note that the semi-thickness H ex rVTr in (|56|). Analogous formulae can 
be found in many monographes (see, e.g., [6], [18], [21] and references therein). 
Thus, in the case of isothermal flow we have 



H oc r 1 " 5 



(57) 



3 Adiabatic flow 



In this section we intend to find an exact steady-state solution to the MHD 
system (Tr5"]) - (122[) provided R e —> oo, i?fe e = 0, and 



0, B r = 0. 



(58) 



Let us remind that the equality v r — implies v z = 0. It could be easily seen 
from the steady-state version of (fT5|) . We take p* = Kpl and, hence, we obtain 
from (JT5J) that 

P = p 1 . (59) 

Then, by analogy with Sec. [5] we obtain: 



B z = B z (r) , v v = v v (r) , B v = B v (r) 



(3 r r 



^_ £ ^ = __ ^ p 



K B, 2 + S3 



(60) 
K97' (61) 



9 o 9$ 

dz F r dz 



(62) 



It is assumed that P (r, z) z=H — and, hence, p (r, z) z=H = 0, where 2H (r) 
denotes the height of disk. Let 



7-1 



^VP = W. 

P 



By virtue of (|55|. we rewrite (R)2"j) to read 

5* 1 9$ 



9z F r -E u 9z 



= 0. 



(63) 



(64) 



We find from flU} 



* + JPC = cw 



(65) 



Since p (r, z) z=H — and, hence, "J (r, 2) z _jj = 0, we obtain from (I65[) that 

1/1 1 



* 



F r £;„ Wr 2 + z 2 Vr 2 + H 2 



Hence 



P 



7-1 



_jF r E u VVr 2 + z 2 Vr 2 + -ff 2 
By virtue of (I5T1) and (j6"3")l . we find 



1 1/(7-1) 



fir r 



-pE u 



^ _ 1. ( E u B l + B l \ _ _P_9£ 
dr dr { fi 2 / F r dr ' 



(66) 



(67) 



(68) 



Eq. ([68]) must be valid under all values of z > 0. Since p (r, z) z=H = 0, we 
obtain from (1681) that 



^ 



a B 2 z + Bi 



r dr \ 2 

By virtue of (J63)) and (|69|) we obtain from (|68|) : 



r 



v 



9* 



1 9$ 9 / $ 



r£ u 9r F r £„ dr dr \ ' F r E„ 
In view of ([70]). we have 

«£ 1 9 



1 9 



F r E u dr ^r^TH 2 ' 



1 



r F r 9r Vr 2 + H 2 ' 

Let, for instance, the motion is Kcplerian, i.e. 



(69) 



(70) 



(71) 



(72) 



then, by virtue of (|TTj) . we obtain that 

H 2 = 0. (73) 

Hence, in view of (joT|) . p = 0. 

We can see from (1TO1) that B v (^ 0) is an even function of z, and, hence, B v 
is unstable (see Sec. Q}. Thus, the stable solution is: B v = 0, and, by virtue 
of (|69p , we find that B z — const. Assuming that the flow is electrically neutral 
|21"1. we find from (l24l) that the flow is the vortex, i.e. 



C, 



i\ 



<p 



<p 



r 



(74) 



Let, for the sake of simplicity, 

C% = ^ (75) 

then, by virtue of (|7Tj) . we obtain that 



H = r\/4r 2 - 1 =>■ JT « 2r 2 . (76) 

Thus, in the case of adiabatic flow, we find that 

# oc r 2 (77) 

4 Perfect gas. Pre- assigned midplane tempera- 
ture 

As it can be seen from Sec. [H the vertically isothermal disk will, in fact, be 
totally isothermal under the assumption that the radial components, v r and B r , 
of velocity and magnetic field, respectively, are negligible. Furthermore, the disk 
cannot be thin in terms of, e.g., [6]. Adiabatic disks (see Sec. EJ, in contrast to 
vertically isothermal ones, are more trustworthy. However, thin disk accretion 
must be highly nonadiabatic, as emphasized in [3D]. Because of this, we will 
consider an intermediate case that is between the above-listed two axisymmetric 
flows. 

In this section we will, mainly, deal with thin disks, i.e. it is assumed that 
E u -C 1. We take p* — p^RT* and, hence, we obtain from (fTTj) that 

P = pT. (78) 

It is significant that the temperature at the midplane, i.e. 

T = T\ z=0 = To (r) , (79) 

is assumed to be a preassigned function of r. 
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4.1 Non-magnetic disk 

It is assumed that 

B r = B z = B v =0. (80) 

Let v z = v r = 0, and let R e — > oo, due = 0. In such a case the steady-state 
version of the system (TT5l) - (f23|) takes the following form. 

— =E U ^ + ^^, P = P T, E u ,F r = const, (81) 

r or b r or 

jMr + f7T = °. s^-7^— Ui^ + o^), (82) 

az F r oz -y/r 2 + z 2 r 2r J v ' 

Since the values p, v v , and P are even functions of z, we can write for small 
values of z: 

p = p + p 2 z 2 + . . . , T = T Q + T 2 z 2 + . . . , v v = v v q + v v2 z 2 + . . . , (83) 

where the coefficients p , p 2 , To, T 2l v v o, and v v2 depend on r, only. Let us note 
that namely z/r <C 1 and, hence, the power series in (1831) must be represented 
in z/r. However, for simplicity sake, the denominators (r 1 , i = 0,2,...) are, 
manely, included into the coefficients. 
Let 

pU«=0, T\ z=H =i). (84) 

Then, by virtue of (|82|) - (l84l . we obtain the following equalities: 

ToP2+T2Po+ 2E^ = ' (85) 

T + T 2 H 2 » 0, (86) 

p Q + p 2 H 2 « 0. (87) 
By virtue of f8l) ]l -([8T l> . we find: 

H 2 « 4T £„F r r 3 , (88) 



T = T 



S^ + ° ^ = r ° - JKKr (^ + ° (<^"» 4 ) ' < 89 > 
P = * (l - i) + O M , P„ - j^ (;) 2 + O ((,/,)*) . (90) 

By virtue of (f83| and (f8T|) . we find the following differential equation in the 
function p — p (r) . 

P -^ = E u ^ + -^,r>r , (91) 

r or r A r r 
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where Tq = To (r) is a preassigned function. We will use the following notation 
Tq = To (ro). The boundary condition is the following 

Po = Po (ro) ■ (92) 

The degenerate equation corresponding to (|91[) . i.e., the equation obtained 
from (EJTJ) by putting E u = 0, will be fulfilled if, and only if, the motion at the 
midplane is Keplerian: 

v v0 =±-±=. (93) 

\Jrt T 

Notice, any arbitrary function p — p (r) fulfills the degenerate equation and, 
hence, there is no problem to take a suitable p (r) such that (l92l would be valid 
and Po ( r ) ~~ * as r — > oo. 

Let E u > 0, then the solution to (JHTJ, ([S2) is the following. 

« [1(2 1 



Po = ^exp / - in j r ( ^o " — ) dr- (94) 



It must be emphasized that p is not an analytic function of the small parameter 
E u , as it can be easily seen from (|M| . Such a situation is common for singular 
systems. Notice, be the disk Keplerian, we would obtain, by virtue of (|94|) . that 



Po = |t\, P°o,Tg = const. (95) 

T o(r) 

Thus, if po M -> as r -> oo, then, in view of (|95|). T (r) — > oo as r — > oo, and 
vice versa, i.e. if To ( r ) — > 0, then p (>")-> oo as r -> oo. Assuming that Tq = 
Tq = const, we obtain that p = p^ = const too. Hence, the assumption that the 
motion is Keplerian leads to improbable density and temperature distributions 
at the midplane under r > r$. To avoid such unlikely solutions, we have to 
accept that the motion is not Keplerian and, hence 

v lo<^jF> r^oo. (96) 

Hence, in view of (|94| , p — » as E u — >• in contrast to the degenerate solution, 
where p (r) is an arbitrary function. Hence, the solution is not a continuous 
function of the small parameter E u . It is also a characteristic feature of singular 
systems. 

It is well known (see, e.g., [5], [20], [22]) that a gaseous disk formed around a 
central star is in an almost Keplerian rotation with a small inward drift velocity. 
Let us investigate this almost Keplerian rotation. Let v r ^ 0, then, instead of 
(l8T |) -(f82 |) . we should consider the following system. 

19(^v) + ^) =0) (9?) 

r or az 
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r * r ^ + Tz {pVrVz) ~ ~f = "5F (KP) - ^ aF' (98) 

- — r /w^ + tt /w^* + ^^ - 0, 99 

r or oz r 

llrWtsWl-sM-^ doc) 

For the sake of simplicity, we take the characteristic quantity, Z» , for length such 
that 

r = 1. (101) 

Since v r is an even function of z, but w z is an odd one, we can write for small 
values of z: 

V r = VrO + V r 2Z 2 + ..., V z = V zl Z + V z3 Z 3 + . . . , (102) 

where the coefficients depend on r, only. By virtue of (l83l) . (|102l) . we obtain 
from dHII), dMI that 

1 d (rp v r0 ) 



dr 



p v zl = 0, (103) 



-■^- r (Po^Vo^ro) + PqV v qV z i H = 0. (104) 

r or r 

By virtue of (TlU3>(fTU4]) . we obtain: 

Po v ro — ^=0. (105) 

ror 

Thus, we find the following exact solution for the midplane value, v v q, of circular 
velocity: 

v^q = —¥-, C v0 — const. (106) 

For the sake of convenience, let us represent C^o as a function of the Froude 
number, F r . Let the point r = r m be the only point where the motion is 
Keplerian, i.e. 

? = i^^£. (107) 

1 m V ' va r r V r 

Hence, 



v v o = v v \ a = — =, r m = const. (108) 

r\J b r 

Notice, in the case of a non-viscous flow with v r o ^ we have the only solution, 
(|106p . for the midplane circular velocity, v v q. If, however, v r o = 0, then we 
have infinitely many solutions for v v o- The only limitations are the boundary 
conditions. In particular, the function w v o = Vipo ( r ) m ust be such that p (r) — > 
as r — > oo in (|94|) . Thus, if v r o — 0, then we have the ill-posed problem, i.e. 
the problem is not well-posed in the sense of Hadamard [T]. It is important 
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to note that any solution (excluding the vortex) for v v o can not be obtained 
as a limiting case (namely, as ivo ~> 0) of the motion under u r o ^ 0. Hence, 
any solution that does not coincide with the vortex will be unstable, as any 
infinitesimal variation, Sv r o ^ 0, gives rise to a finite response in the gas flow. 
Nevertheless, the power-law model [351 P- 374] 



a 



f 



C v = const, 0.5 < ^ < 1, (109) 



is assumed to be stable (Rayleigh stable [7]) to pure hydrodynamic perturba- 
tions, since it satisfies the so-called Rayleigh stability (necessary and sufficient) 
condition [21 p. 78]: 

fl^^l^) 2 >o, (no) 

where tt denotes the angular velocity, f\ denotes the Rayleigh frequency. 

Let us estimate v V 2 in (JH3J)- We will use the following asymptotic expansion 

oo 

To (r) « J2 a " r ~ n (HI) 

71=0 

in the limit r — > oo. As usually, the series in (jllip may converge or diverge, but 
its partial sums are good approximations to To (r) for large enough r. Assuming 
that To (r) — > as r — ¥ oo, we find that ao = 0. Then, assuming the character- 
istic quantity, T*, for temperature such that a\ = 1, we represent To (r) in the 
following form: 

T (r) = - + O (r- 2 ) , r -> oo. (112) 

Thus, we can use the following approximation: 

T (r)«~, r->oo. (113) 

r 

One can easily see, for instance, that the midplane temperature of the adiabatic 
disks considered in Sec. [3] can be approximated by (J113I) . 
By virtue of ([TT3| . we find from ([88]) that 

Hocr. (114) 

Let v z is a linear function of z, i.e. v z = t> z i,z with u z i = v z \ (r). Since 



Vr 



u z\z=H ~ Q r Ur \z=H ■> 

we find, by virtue of (fTT3|) and ([88]), that 



(115) 



«2l-« ~ ^— «r0 =► V zl « . (116) 

or r 
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By virtue of ([TU21 and ([55]). we obtain from ([571), ® that 



^0-^ + 2^1^2 + -^ =0. (117) 

ar r 



Thus, by virtue of (|116p . we find 



const ,„ „. 

v v2 « — «-. (H8) 



By virtue of ([83 ]) . ([88 ]) . (fTT3l) . and ([TT8l) . we find that 



v v \ z=H ™—(l-C v2 E u Fr), C v2 = const > 0. (119) 

Since -B u 4C 1, it can be assumed that 

V V « Se°, < z < ff. (120) 

It is easy to see from (|120p that the Rayleigh frequency /£ = and, hence, 
the vortex is Rayleigh unstable. Let us consider the problem of vortex stability 
from another point of view. 

4.1.1 Stability 

Let <T denote the value of a dependent variable <j for the case of a steady-state 
solution. It is assumed that we consider the motion close to the solution to the 
system ffHTj) - (f82]) . In view of (|108[) and (|113[) . the solution can be written in the 
form: 

v z = v r = 0, v v0 = v v \ z=Q = — -=p=, r m = const, (121) 



rVF r 

| = rexp /_L(^_J_),, (122) 

l 

In such a case the non-linear system in perturbations ? (<j = <T + ?) for 
T5)) - (rT5)) will be, under _R e — > oo, as follows. 

a 9{p + p ) , l dr{p Q +-p )v ra 

di r d~r ^ Po Pa > zl ' ^ ^ 



c 9(Po+ft)"rtl9 ,_ .~x~ 2| ,- L ~x~~ (p + Po) ( v vo + v v0 ) 2 
ot r or r 

-E u §- r (p ^ )(T + f )-^, (124) 

Sh d(Po + ft) ft* + V>) + Id 

ot r or 
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(Po + Po) ( v ^o + v v q) Vzi H = 0, (125) 

c, g (Po + Po) ^1 19 /_ . ~ s ~ ~ . ,_ . ~ x ~2 

at + r* r ^ Po Po ^ UzlWr0 ^ Po Po ^ Vzl = 



2E° r 

oz 



(Po + p ) (r a + r 2 ) + ( P2 + p 2 ) (t + f )] - ^±|^ (126) 

Our main purpose is to reduce ()123|) - (|126p to an ODE system for a subse- 
quent stability investigation. We will consider a function as admissible if it will 
be a singular function |13j . That is, this function is continous on 1 < r < oo 
and the derivative over the coordinate r exists and is zero almost everywhere. 
It can be also assumed that the function is strictly monotone decreasing, e.g. 
Vr\ r - j.oo ~* 0- Cantor staircase [13] and Lebesgue singular function [11], [12] 
are well known examples of such functions. Singular functions occur in physics, 
dynamical systems, etc. (see e.g. the references in [TT]). 

We restrict ourself to the case of linear system. Assuming the non-linear 
terms in f|123|) - (1126[) as negligible, we get: 

dp 1 drp v r0 _ 

h ~dt r dr + PoVzi-U, (127) 

_ dv r0 2-p v v0 v V Q + p vl 
ShPo ^l r " 

-E u ^-^ T + p T )--^-, (128) 

^f = 0, (129) 

ShPo ^f - -2E U (p f 2 + p Q T 2 + p 2 f + p 2 T ) - A . (130) 

By virtue of dHHJ) - dHOJ) , we obtain from (1X27]) - {150]) : 

d'Po ~ , / 1Q1 -> 

— — = ai2V r0 +ai3V z i, (131) 

-rrr = a 2 {p Q + b 2 , (132) 

dvzi 
dt 
where 



h, (133) 



drPo _ Po C lo j? dT Q 1 

12 = ~7?I7P 13 ~ ""SV ' 21 " r3S*,75 ~ u 9,fl flr ~ 75 r 2 S,F ' ( j 

I OhUI Oh I OhPQ DhP(jUI PqI Dh r r 

2C v0 ~ „ <9p ~ 2T 2 , 1Q ,i 

02 = -oTT^O - E u _ T , 63 = -£„-£-. (135) 

H5/, S h p dr S h 
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Hence 

d 2 ~ Po 



ai 2 a 2 ip = ai 2 b 2 + 01363. (136) 



dt 2 

Obviously, Eq. (|136[) and, hence, the system (|13ip - (|133[) will be unstable if 
ai 2 a 2 i > 0. 

By virtue of (fT2"2l and ([TT5|) . we find 

E| = r .-,„p(f,--f), ,- a ^>I,^<l. (137) 

In such a case we obtain that 012021 > if 

[ FrC^o C 2 F r \ 

r > max — , . (138) 

yi-E u F r , (2E u F r -l)J v ' 

Since the disk is thin, i.e. E u <C 1, and using (|107[) , we find the following 
sufficient condition of linear instability 

r > j-^f-p , (139) 

i — rj u r r 

where r m is the only point at which the motion is Keplerian. 

Thus, the vortex will be rather unstable, at least, on the disk's periphery, 
provided R e — > 00. 

4.2 Magnetic accretion disk 

It is assumed that all dependent variables are analytic functions of z in a vicinity 
of the midplane. Then, in view of the symmetry, we may write: 

B v = B v iz + . . . , B z = B z0 + B z2 z + . . . , B r = B T \z + B^z + . . . , 

P = Po + P2 z2 + • • • ) T = T + T 2 z 2 + . . . , v v = v va + v v2 z 2 + . . . , 
IV = VrO + v r2 z 2 + ..., v z = v zl z + v z3 z 3 + .... 

$ = - J—, = $0 + ®2Z 2 + ... = -- + ^z 2 + ... . (140) 

Vr 2 + z 2 r 2r A 

It is also assumed in this section that B M 7^ 0. 

A steady-state magnetic accretion disk with v r o 7^ 0, p 7^ 0, and with 
a negligible dynamic viscosity, i.e. /i = 0, will be our initial concern. The 
magnetic field will be called as almost poloidal if 

B vl = =*• B v = O (z 3 ) , (141) 

and the magnetic field will be called as almost axial if 

B vX = 0, B rl =0^B v = O (z 3 ) , B r = O (z 3 ) . (142) 
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Let us prove that the circular velocity at the midplane will be the vortex, 
i.e. 

v v0 = v v \ z=0 = —*—, C v0 = const, (143) 

if and only if the magnetic field will be almost poloidal. Actually, since p, = 0, 
the steady-state version of Eqs. (fl"5j) . (fT7|) is the following. 



1 d (rpv r ) d (pv z ) 





r dr dz 




1 d ( E " R R V 

r dr~ T \ P v r ~ ~J B * Br ) + 


d 


/ E u \ pv v v r E u B (p B r 


dz 


\pv v v, p B v B z j | r p r 



(144) 



(145) 



In view of (|144p . (|145p . and (|140p . we obtain 

l<9(rp w r0 ) 



<9r 



Po^i = 0, (146) 



Sr r (Po v ro)-\ ^r ^-B vl B zQ +p v v0 v zl + - =0. (147) 

r or r or p f 

By virtue of (|146l) . we obtain from (|147|) : 

Po v r0^r ^-B vl B z0 + ^ = 0. (148) 

or P t 

Let the magnetic field will be almost poloidal, i.e. B v \ = 0. Then, in view of 
(fT48|) . we have 

Po^o^=0. (149) 

ror 

Equality (1T4T))) proves ([Lt3"l) . 

Let (J143P be valid. In view of (J148I) . we have 

Po«rO — ^ ^-B v iB z0 = 0. (150) 

rar p 

By virtue of (fT43|i . we find from (|T50| : 

S v i^o - 0. (151) 

Equality ([T5H proves ([CTj) . 

Let us now consider a steady-state magnetic disk with v r = (and, hence, 
v z = 0) and with a negligible dynamic viscosity, i.e. u = 0. Let us remind 
that such assumptions for the case of non-magnetic disk lead to the ill-posed 
problem in the sense of Hadamard (see Sec. I4.ip . In particular, the solution to 
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the mathematical model (Sec. I4.1[) is not unique. We intend to find a steady- 
state solution to the MHD system (ll~5l) - (j22l) provided R e — > oo, -dke = 0. In 
such a case, the MHD system is reduced to the following. 

E U B% pvl_ d ( E U B*\ p3$ 

-\HT B » B <)-t{T B » B )-r-^=°' «-» 

d{v v B r ) , d(v v B z ) =Q) (155) 



<9r <9z 

1 drR. <9£ z 



r dr dz 

Eq. ([24"]) is reduced to the following 



= 0. (156) 



dT{v f z) ~ d{V f r) = 0- (157) 

ror oz 

Let B v be an odd function of z and let B r = O (z 5 ) . Then we find from (|156|) 
that 

B z0 = B z0 (r) , B, 2 - B z4 = 0. (158) 

By virtue of ([T551) . we find from (|T55l) and (fT55f that 

u^o = %o {r) , v v2 = Vipi = 0, B v i = B v3 = 0. (159) 

Taking into account that B r = O (z 5 ) and using (I158p . we obtain from Eq. 
(|154j) that (|85|) is valid for the case of the magnetic disk as well. Hence, using 
dHSl), (J901), and (fi~58 l) - ([T59l) . we obtain from Eq. (fT52f that 

^ = E,T% + Ej-L + ^ + JL?l + 0( r J). (160) 

r or or p or b r or v ' 

Eq. (|160p must be valid under all values of z. Assuming z = H, we obtain, in 
view of (|84|). the following equation in B z q: 



E u dBl a 
/3 dr 



O 



((z/r) 4 ) = 0. (161) 
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Then we find: 

B z q w const. (162) 

Since B r = O (z 5 ), we obtain from (fT5?[) . by virtue of l[T5g|l - ([T5g|) . that 

rar 

We find, in view of |[IB2|) . ([IBS]) , that 

const „„,. 

%0 a . 164 

r 

Analogously, if B r = (and w r = 0), then we find from (|152[) - (|156[l that 
B v = and B z = const. Thus, in view of (|157l) . the vortex 

const 

v v = 165 

r 

will be the only solution for the circular velocity, and, hence, the motion will 
be Rayleigh unstable. Let us note, the density and temperature distributions 
can be easily found from Eqs. (|152|) . (|X54|) . which can be written, in view of the 
foregoing, as the following: 






(E uPo T Q ) + -£^, (166) 



where the circular velocity, v v , is calculated from Eq. (|165|) . Notice, Eqs. (|166|l . 
(fWj) coincide with Eqs. dHU), O, and, hence, we obtain ([55]), (1501) . and ([Ml) . 
Notice, we have assumed v r Q ^ in the above-proven assertion that the 
midplane circular velocity will be the vortex if and only if the magnetic field 
will be almost poloidal. The following counter-example demonstrates that the 
condition v r o ^ is essential. We consider the case when v r = and B z = 
B z o + O (z 4 ) with B z o — const. In view of (|156p . we have 

C 

B r \ = , C r \ = const. (168) 

r 

Then, by virtue of f|168[) . we obtain from (|157[) the power-law model: 

const _ C r i 

v v q = , x =1 = const, (169) 

r* B z0 

which is the only solution for the midplane circular velocity in the frame of our 
assumptions. However, in view of (|153[) . B v \ = 0, i.e. the magnetic field is 
almost poloidal. 



20 



4.2.1 Viscous disk 



Let us, first, demonstrate that the vortex, (|143p . will also be the midplane 
circular velocity if the dynamic viscosity fi = const ^ 0, provided that the 
magnetic field will be almost axial. Actually, if B v = O (z 3 ) and B r = O (z 3 ), 
then, in view of (|140[) , we obtain from the steady-state version of Eq. (fT7|) that 



drv, 



Po v rO- 



'</>o 



rdr 



J_d_ 

BZlh- 



fir 



#(-) 
or \ r / 



R e R e dr\r) 



In view of ([221) and (fT4T)|) . we have 

1 9rB r i 
r dr 

Hence 



2B z2 = 0. 



5,2 = 
By virtue of ([MOf . we find from d20jt that 

9 (iVo-B^i -^o-Sri) 



2 (v zl B vl - v v2 B z0 - v v0 B z2 ) = 



Hence, in view of (|172[) . we obtain: 

v V 2B z0 = =*• v v2 = 0. 
By virtue of ([T70]) and (JT741 . we find that 



Pa v rQ- 



drv, 



Vo 



r(9? 






d /v, 



ch 



(*?) 



,V> 



(170) 

(171) 
(172) 

(173) 

(174) 

(175) 



Obviously, the vortex, (I143|) . fulfills Eq. (I175|) . To estimate other solutions, v v q, 
to Eq. (fT75]> . we will use (fTTo) . i.e. 



^ = ^ + 0((H/r) 2 ). 
By virtue of (|176p . we obtain from (J103I) that 



r<9r 



Hence, we find with accuracy O ( (H/r) ) that 



Qo 



rp v r Q = , Qo = const > 0. 



(176) 



(177) 



(178) 



Since /i = const ^ 0, we assume that jj, = 1, without loss in generality. Then, 
by virtue of (|178[) , we rewrite (|175|) to read 



a r dr 2 £lo d ( <9f2o 



d) 



dr 



2J7q , a r — QoR c = const, 



(179) 
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where fio denotes the angular velocity at z = 0. We will investigate (1170)) in 
the vicinity of r — 1. Let x = lnr, then we have instead of (|1T9[) 

-—f- + [2 + <v exp (-a;)] -^- + 2a r exp (-x) O = 0. (180) 

In the vicinity of r = 1, i.e. in the vicinity of x = 0, we obtain, in view of (|180[) . 
that 

«o«^ + % Ci,C 2 = const. (181) 

Thus, if a r > 2, i.e. if i? e will be sufficiently big, then the motion will be 
Rayleigh unstable, and, hence, the motion can be stable under sufficiently small 
Reynolds number. In all likelihood, it will be a turbulent motion. 

As indicated above, the midplane circular velocity can be the vortex if the 
dynamic viscosity fi = const. Let us now consider the case when the midplane 
circular velocity is not the vortex, but the power-law model: 



a 



'J<pO 



•p 



const t^ 1, Cp = const. (182) 



The power- law model, (fl8"2]) . can fulfill Eq. l[T?5|l if p, = fi(r). The following 
inequality must be valid to fulfil (IM1) . 

x>0.5. (183) 

It is also assumed that the magnetic field will be almost axial, i.e. B v = O (z 3 ) 
and B r — O (z 3 ) . To estimate fx = \x (r) at z = we assume v r o = 0. We obtain 
from the steady-state version of Eq. (|17|) that 






dr V r I 



R e dr\ r ) 



0. (184) 



Then, by virtue of ([152^) . we find from ([TM)l that 



C 

/i = 1 _ A '^ , C M = const. (185) 

Assuming the characteristic quantity, u rl for viscosity such that C M = 1 and 
using (|113|) . we represent the viscosity, /z, in the following form: 

H = T 1 -*. (186) 

The power law (e.g., [4], [16], [23]) for the laminar viscosity, fi h of dilute 
gases can be written as the following 

Hi = T\ (187) 

where typically 6 = 0.76 |4J. It is, also, assumed that 6 = 8/9 if 90 < T < 300 
°K, and 9 = 0.75 if 250 < T < 600 °K. If we assume that the flow in question 
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is laminar, then /i = /i ; , and, hence, x. = 1 — 9. Since we consider the flow at 
the periphery of the disk, i.e. under a low temperature, we find that 

x < 0.25. (188) 

The inequality (J188I) contradicts to (I183p . and, hence, the flow in question is 
rather turbulent. 

The Euler number characterizes "losses" in a flow [10] . and it is higher in a 
turbulent flow than that in the laminar regime [5] . Let us now estimate the value 
of k in (j!82[) by assuming that the turbulent gas tends to flow with minimal 
losses, i.e. to have the Euler number as small as possible. Using Prandtl and 
Kolmogorov suggestion [TJ p. 230] that the turbulent viscosity, n 1 is proportional 
to the square root of the kinetic energy of turbulence, k, we evaluate u = p\ z=0 
as 

fi Q = C K L K p K - 5 , C K ,L K — const. (189) 

By virtue of (|TT3"|) . (|TUTj) . and (TIE?)) , we rewrite flM]) to read: 



,fl 



/ ? o = -^T ex PC 



1-1 



£;„ (2x-l .E„.iv 



Let the kinetic energy of turbulence k = const. In such a case, we obtain from 
CPU) that 



^o = — expC 



Cp = p° C K L K T^- 5 = consi. (191) 



Equating (|185[) and (|19ip at r = 1, we find that C p = C M . Let us now assume 
that /x of (|185|) and /i of ()191|) coincide each other in the vicinity of r — 1, i.e. 
at r = 1 + e (e <C 1), with accuracy O (e 2 )- In such a case we obtain that 

£„oc— !— . (192) 

Z — H 

As we can see from ()192|1 and (|183j) . £"„ — > min if x — > 0.5. Thus, in the frame 
of our assumptions, we find that the turbulent flow tends to be Keplerian. 

5 Concluding remarks 

On the basis of the exact solution to the MHD system in Sec. [5]we may conclude 
that vertically isothermal disks will, in fact, be totally isothermal under the 
assumption that the radial components, v r and B r , of velocity and magnetic 
field, respectively, are negligible. Furthermore, the disks cannot be considered 
as thin in terms of, e.g., [B] even if the Euler number < E u 4C 1. 

The exact solution to the MHD system in Sec. [3] corroborates the view [5D] 
that thin disk accretion must be highly nonadiabatic. Despite of the fact that 
adiabatic disks (see Sec. [3j) are more trustworthy than isothermal ones, we find 

ocr na. 
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The exact solutions to the MHD systems in Sec. [5] and Sec. |3] prove that 
the vortex will be the only solution for the circular velocity provided that the 
flow is charge-neutral. Let us note that the exact solutions are found under the 
assumption that v r — and B r = 0. 

The approach developed in Sec. 2] for the modeling of thin accretion disks 
turns out to be efficient. In the case of non- magnetic disk, this approach enables 
to obtain, with ease: the solution for the steady-state non-viscous disk with a 
good accuracy, to find the non-dimensional semi-thickness H ex r, to prove that 
all solutions for the midplane circular velocity are unstable provided the disk 
is non- viscous. Using this approach one can prove that the midplane circular 
velocity will be the vortex if and only if the magnetic field will be almost poloidal. 
The approach of Sec. [4] enables one to demonstrate that the pure hydrodynamic 
turbulence in accretion disks is possible, and to demonstrate that the turbulent 
flow tends to be Keplerian. 
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